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ABSTRACT 



This thesis examines the problem of predicting distortion 
induced stall' in an axial flow compressor. Three instantane- 
ous distortion indices are used on the data presented. The 
data consist of eleven stall events run at 2.5 Mach by NASA 
Lewis on a J85-GE-13 turbojet engine. The three indices 
used include two indices based on vorticity and developed at 
the Naval Postgraduate School. The third index was developed 
by Pratt and Whitney and is based on total pressure distri- 
butions. The three indices are compared by visual means us- 
ing computer plots and are also compared by using a simple 
statistical approach. The computer program used was first 
developed by Shoemaker in his Master's thesis. The program 
is expanded in this thesis to compute and compare the three 
indices . 
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I. INTRODUCTION 



This paper addresses the problem of distortion-induced 
or drift stall in an axial flow compressor. The problem is 
a current one in industry and has not been wholly resolved. 
The present method of dealing with the problem of distortion 
is to operate axial flow compressors with a sufficient stall 
margin to avoid the problem altogether, but this is an expen- 
sive procedure. Best engine efficiencies are obtained by 
operating close to the stall or surge line. By providing a 
large stall margin, we are in effect sacrificing performance 
via reduced efficiencies. The problem is especially germane 
to military aircraft. The normal operating conditions of 
military aircraft include hard maneuvering, high "g” flight, 
steam ingestion from aircraft carrier launches, and exhaust 
gas ingestion from forward firing weapons. These operating 
regimes cause a pressure, temperature, and gaseous molecular 
weight distortion of the flow entering the compressor section 
of the engine. 

In order to deal with distortion, it is first necessary 
to find or define an adequate measure of it. A distortion 
index based on steady-state or long-time averaged values of 
distortion has been attempted. This method has met with some 
degree of success, but it is apparent that peak instantane- 
ous distortions of sufficient magnitude can induce stall when 
the time averaged distoi'tion factor remains relatively low. 
The current approach and the approach of this paper is to 
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examine distortion as an unsteady phenomenon and to measure 
distortion over time intervals comparable with the character- 
istic times of an axial compressor. 

A computer program was first developed by Ens. Shoemaker 
[ref 3] utilizing an existing interpolation scheme developed 
at the Lewis Research Center of NASA. The program calculates 
the total pressure distribution just prior to the compressor 
section of a J85-GE-13 turbojet engine as a function of time. 
The program also utilizes the theory of Fuhs, Farmer, and 
Iverson [ref 1, ref 2] to compute the circumferential and 
radial components of vorticity across the face of the com- 
pressor section as a function of time. The computer program 
was extended for this paper to calculate and statistically 
compare two existing distortion indices plus an additional 
index developed by Downs and studied in this thesis [ref 4]. 
This procedure was applied to eleven stall events recorded 
and analyzed by the NASA Lewis Research Center [ref 5]. 
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II. NATURE OF THE PROBLEM 



A. COMPRESSOR SURGE 

The line marked "surge" in figure 1 denotes the locus of 
unstable operation points for an axial compressor. The 
cause of compressor surge may be visualized in the following 
manner. At any given reference RPM (N//F) and blade radius, 
a rotor speed U is specified. Consider the velocity diagram 
in figure 2. A reduction in referred mass flow rate will 
cause a reduction in the axial component of velocity. Since 
the rotor speed is fixed, an increase in the relative flow 
angle 3 will occur. Present compressor design incorporates 
a fixed rotor geometry and thus the relative flow angle may 
only increase to a maximum value which is dependent on blade 
design . 

An insight into the unstable nature of the surge line 
may be made by comparing the characteristics of the compres- 
sor to the characteristics of a throttle. A throttle is set 
to produce a certain back or receiver pressure, given an in- 
let pressure. We may visualize a system consisting of a 
compressor producing a high inlet pressure followed by a 
throttle which dumps the high compressor air to ambient con- 
ditions (fig 3). A plot of the characteristics of the hypo- 
thetical throttle and compressor are included in figure 4. 
The compressor is represented at a fixed reference RPM. The 
dashed lines represent fixed throttle settings. Selecting a 
particular throttle setting determines two possible 
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Compressor- Throttle System 



FIGURE 3 




FIGURE 4 
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compressor points. Consider as a throttle setting line B, 
with possible operating points P and Q. First consider point 

P. If the mass flow rate is perturbed to a higher mass flow 
rate, the throttle setting will demand a higher pressure 
ratio, whereas the compressor produces a lower pressure ra- 
tio, In this case the combination of compressor and throttle 
will return to point P. Conversely a lower flow rate will 
cause the throttle to demand a lower pressure ratio whereas 
the compressor will produce a higher ratio. Again the sys- 
tem will be forced to return to point P. Consider now point 

Q. A slight perturbation mass flow will produce problems. 

The system is free in this case to follow any perturbations 
and the system is unstable. Neutral stability is represent- 
ed by point R where the slope of the throttle curve matches 
the slope of the compressor curve. 

B. CLEAN STALL VERSUS DISTORTION INDUCED STALL 

The surge line in figure 1 denotes the conditions at 
which stall will occur if the inlet velocity profile is uni- 
form, producing a uniform total pressure distribution. In 
actuality, the velocity profile is not uniform and pressure 
perturbations may produce a local condition over a small re- 
gion of the compressor face which will cause stall over that 
region. This stalled region in turn produces a condition 
known as rotating stall. In figure 5, if the conditions are 
such that blade 2 becomes stalled, the effective flow area 
between blades 2 and 3 is reduced. This reduction in flow 
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Rotating Stall 



FIGURE 5 



area causes a flow diversion which tends to overload and 
stall blade 3. Continuation of this process will produce 
rotating stall. It is noteworthy that rotating stall propa- 
gates in the opposite direction from rotor velocity at about 
half rotor speed. Rotating stall may develop into surge or 
fully developed stall. 

It is apparent that the instantaneous velocity distribu- 
tions across the compressor face may cause stall when the 
steady-state conditions are below the surge line. This ve- 
locity distribution in turn produces a total pressure distri 
bution which is the form of the data examined in this thesis 
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III. DISTORTION INDICES 



A. VORTICITY VIA CROCCO'S THEOREM 



Non-dimensional vorticity components in both the circum- 
ferential and radial directions of a cylindrical coordinate 
system may be calculated using the theory of Fuhs, Farmer and 
Iverson [ref 1, ref 2]. Starting v/ith Crocco's Theorem and 
using an order of magnitude analysis, it was possible to re- 
duce the expression for these parameters to: 

, _ -1 ap' M 3 t' . , _ 1 ap' M aT' 

“e yP'U' 3r' 2 ar' ’ “r yP'r’U' ae 2r' 36 
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For a detailed explanation, consult reference 1. The calcu- 
lations are examined more briefly in reference 2. 

Analytical studies and data analysis have been performed 
in an attempt to relate the effects of total temperature and 
total pressure distortions on compressor performance [ref 7] 
The vorticity approach is well suited to this problem, as 
both total temperature and total pressure may be taken into 
account in the calculations of vorticity. It would be bene- 
ficial to examine data which consist of both temperature and 
pressure patterns across the compressor face. Additional 
work in this area would be helpful. The data examined in 
this thesis, however, did not present a significant tempera- 
ture gradient and the expressions for vorticity could be 
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reduced to: 




= -1 9P 



1 ap' 
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The task of this thesis was to determine whether an index 
based on vorticity would provide a meaningful prediction of 
compressor stall. Two such indices based on vorticity are 
examined in this thesis. The first index was proposed by 
Shoemaker [ref 3] and is based on circumferential vorticity. 

A second index proposed by Downs [ref 4] and based on both 
radial vorticity and total pressure is examined. In order 
to provide a means of comparison, a third index based on to- 
tal pressure and taken from Zonars [ref 8] was also examined. 

B. SHOEMAKER INDEX 

James E. Shoemaker approached the problem of finding a 
vorticity based index by first printing out maps of pressure, 
circumferential vorticity and radial vorticity. In order to 
accomplish this goal he utilized a set of subroutines devel- 
oped at V/right Patterson Air Force Base, and altered them to 
print out the maps at times near stall and at times far from 
stall, to provide a comparison. A copy of his computer pro- 
gram is included at the end of his Master's thesis [ref 3]. 
Using this procedure. Shoemaker found a vorticity pattern 
that frequently occurred within twenty-five milliseconds of 
stall, and was not common at other times. The pattern was 
characterized by a ring of large negative circumferential 
vorticity at the midspan of the compressor face, and a large 
positive ring at the tip. He developed an index based on 
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this pattern which will be referred to as the Shoemaker in- 
dex. He divided the compressor face into three concentric 
regions. Region one at r' less than 0.5 was ignored. Region 
II was established as that area between radii of 0.55 and 
0.70. Region III was that area between radii of 0.8 and 
0.95. The Shoemaker index then summed over Region II but set 
any positive circumferential vorticity equal to zero. Like- 
wise Region III was summed where any negative vorticity was 
set equal to zero. The instantaneous Shoemaker index was 
set equal to the sum of the absolute value of the two summa- 
tions. In notation form, 

Shoemaker Index (T) = E ^ E [Wg] 

where = positive values of circumferential vorticity in 

b 

Region III 

to„ = negative values of circumferential vorticity in 
b 

Region II 

In order to take into account the added effect of the 
duration of the pattern, Shoemaker set his final index equal 
to the product of three consecutive instantaneous indices de- 
scribed above. This latter procedure was not used in the 
present thesis. The Shoemaker index as used herein thus re- 
fers to the instantaneous index and is an indication of dis- 
tortion in the radial direction. 

The Shoemaker index has one real advantage when the 
transport time from the measuring probes to the compressor 
stage (where stall originates) is considered. The pattern 
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develops at a time before stall which is larger than the 
transport time from the measuring probes to the stall region. 
This could be very useful if variable stators or the fuel 
control were available to avert stall. 

C. DOWNS INDEX 

A second index which was proposed by Downs [ref 4] was 
formulated and is explored in this thesis. This index incor- 
porates both total pressure and radial vorticity. The basic 
idea was to consider the contribution of both total pressure 
and vorticity to the relative flow angle 3- Increasing 3 
would cause the rotor blade incidence angle to increase and 
would tend to load the blade. A perturbation plot of 3 ver- 
sus 0, where the deviation of 3 from the reference line takes 
into account pressure and vorticity effects, is shown in fig- 
ure 6. The corresponding plot of total pressure and is 
also shown. The reference line represents uniform flow. It 
can be seen that although is dependent upon the pressure 
gradient, the combined effect of pressure and on 3 is de- 
pendent on both the gradient and relative magnitude of pres- 
sure. The Downs index was formulated to search the compressor 
face for the largest value of the combined effect of pressure 
plus over any 60 degree sector. In notation form, 

P -P 

Downs Index (t) = Z [— ^ — ^ + to'] 

P,p ^ max 60 

where P,j, = average total pressure over entire face of 
compressor at time t 
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A = possible 
B = corresponding 
C =- ^^-^ -plus u' effect on B 
D effect on B 

0 measured in direction of rotor rotation 

Effect of pressure and radial vorticity on (3 

FIGURE 6 
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The Downs index is primarily an indication of distortion in 
the circumferential direction. For a more complete analysis, 
consult Downs [ref 4]. 



D. PRATT AND WHITNEY K„,, DISTORTION INDEX 

DM 

In order to provide a means of comparison between vortic- 
ity based indices and an index calculated directly from total 
pressures, the Pratt and Whitney index was selected from 

Zonars [ref 8]. This index was selected primarily because 
of case four, which was one of the eleven cases investigated. 
Distortion indices are plotted for all eleven cases in Appen- 
dix A. Neither of the vorticity indices had satisfactorily 
predicted stall in case four; from examining the pressure 
maps, Kj^j, appeared to be an appropriate index. The index 

was applied to case four and to all other cases. This index 
was developed by Pratt and Whitney for the TF-30 installed in 
the F-lllA aircraft and is defined as follows: 






where 



P - P 

1 ^ r Tmax Tmin, ^ 

2 E [ p- ]. e C. 

Tave 

E C. 

1 



X 100 



C = ratio of compressor inlet radius to ring radius 
i = number of ring 

9~ = largest continuous arc of the ring over which 
the total pressure is below the ring average 
pressure 

^Tmax ~ maximum total pressure 



^Tave ~ average total pressure 
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^Tmin ~ minimum total pressure 

It can be seen that the K^,, index is an indication of dis- 

DM 

tortion in the circumferential direction and puts more empha- 
sis on the region near the hub. It would be expected that 
the index and the Downs index would be comparable since 

both search for circumferential distortion. This was in 
fact the case and will be discussed in a later section. 
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IV. EXPLANATION OF DATA 



The NASA Lewis Research Center conducted an experimental 
wind-tunnel investigation to determine the effects of time- 
varying total-pressure distortions produced in a supersonic 
inlet on a J85-GE-13 turbojet engine [ref 5]. Selected data 
taken during this investigation were digitized on magnetic 
tape; the latter is the source of data used in this thesis. 
The following explanation of data is largely a brief summary 
of the explanation presented by Burstadt and Calogeras [ref 
5] and is given here to provide an overall understanding of 
the data analysis. 

The inlet was an axisymmetric mixed-compression type. 

The compressor was an eight-stage axial flow compressor with 
interstage bleed and variable inlet guide vanes. An after- 
burner section was present but was not used during the test 
procedure. The engine bleed system was modified to make the 
compressor operate with less than normal stall margin. The 
following excerpt is quoted from reference 5: 

The test procedure was to establish an engine 
operating condition in the vicinity of the com- 
pressor stall point. The FM tape recorder was 
started and a steady-state data scan was made. 

Then, while the FM tape recorder was still run- 
ning, the compressor pressure ratio was increased 
slightly and the throttle adjusted to maintain a 
constant rotor speed. If stall did not occur, 
this procedure was repeated until the stall event 
was recorded on the FM tape. In some cases the 
stall occurred after the engine had stabilized 
at an equilibrium operating condition. These 
were called 'drift' stalls. 
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All eleven of the stalls examined in this thesis were of 



the "drift” type. NASA Lewis considered only seven of the 
eleven cases in their analysis. Cases 6, 8, 9, and 10 were 
excluded on the grounds either that their causes were attri- 
butable to steady-state distortion level.s or that they showed 
very little loss in stall margin. Table I shows steady-state 
operating conditions for all eleven cases presented in this 
thesis. The NASA Lewis analysis identifies each case number 
as a steady-state reading number. These NASA identifying 
numbers are included in Table I for comparison [ref 5] . 

Steady-state and dynamic pressure instrumentation at the 
compressor face is shown in figure 7. Only the fluctuating 
component of pressure was recorded from each of the dynamic 
probes. Absolute pressures were later obtained by adding 
the fluctuating component to the steady-state value. 

Frequency response of the dynamic probes was flat to 
about 2000 hertz. The fluctuating signal from each pressure 
transducer was recorded on FM multiplexed tape at 152.4 cen- 
timeters per second. The data in a 225 millisecond period 
just prior to stall were digitized at an effective rate of 
8000 points per second per channel. These data when processed 
by the computer program listed in the back of this thesis 
provided the distribution of total pressure across the com- 
pressor face with a time spacing of 0.125 milliseconds between 
distributions. The distributions of total pressure were fur- 
ther processed to provide distributions of both circumferen- 
tial and radial vorticity as a function of time. In order to 
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Compressor Pace 
FIGURE 7 

« 

reduce the amount of calculating involved, only one out of 
five readings was actually processed so that there was five 
times 0.125 or 0.625 milliseconds of time between readings 
actually processed. 

The detex'minat ion of stall time was made by NASA Lewis. 
The signals from the compressor static pressure transducers 
were put on strip charts and examined by NASA. For all 
cases, the row of instrumentation at the 290 degree location 
showed the earliest indication of stall. Since the origin 
of stall in a compressor stage causes an abrupt loss in the 
airflow pumping of that stage, a compression wave is 
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propagated upstream. This compression wave is referred to 
as a hammershock. The following excei'pt is quoted from ref- 
erence 5: 

Because there were only two axial rows of in- 
strumentation inside the compressor, it was not 
possible to determine exactly when stall originated 
inside the compressor. Although the first indica- 
tion was seen at the 290 degree location, a stall 
zone could originate just past the 45 degree loca- 
tion but not be detected until it rotated past the 
instrumentation at 290 degrees. Since the stall 
zone rotated at about half rotor speed (8 msec/rev) 
there could be as much as a 6-millisecond delay be- 
fore it was detected at the 290 degree location. 

A one-millisecond period for transport time between sensors 

and stall stage, plus time for blade response, was predicted 

by NASA. This produces a six-plus-one or seven-millisecond 

possible time lag between actual initialization of stall and 

the stall time listed in this thesis. 
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V. THREE INDICES APPLIED TO DATA 



A. EXAMINATION PROCEDURE 

The computer program first developed by Shoemaker [ref 
3] was adapted to compute each index as a function of time 
for the entire time interval provided by NASA Lewis. This 
time interval varied slightly but. was generally about 215 
milliseconds of data prior to first indication of stall. 

The DRAW subroutine of the Naval Postgraduate School was used 
to plot each index. These plots are included in Appendix A 
and are plots of the instantaneous version of each index. A 
numerical printout was also provided by the program. A sam- 
ple numerical printout for the Shoemaker index case one is 
included at the end of this thesis. Table II provides an 
explanation of the terms for the numerical printout. 

In addition to the computer plots, which may be used to 
determine if the given index is an accurate predictor of 
stall, a simple statistical approach was used. After the 
first fifty values of the index were computed, the mean and 
standard deviation for those fifty values were obtained. 

Then, as the next fifty values of the index were computed, 
the following test was made. If the average of the three 
previous values of the index at any time exceeded a certain 
criterion value, the computer program printed out a stall 
warning with the time in milliseconds of occurrence. The 
criterion level was set in the following manner: 

Test = mean + 2.5 x standard-deviation. 
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TABLE II 



EXPLANATION OF PROGRAM TERMS 



N 

FACT(N) 

DFACT(N) 

TIME 

PTAVE 

MEAN 



STD- DEV 



= index number (a function of time). 

= value of instantaneous index at time ''N'' . 

= average of three previous FACT(N) values. 

= "N" time in milliseconds. Starts at zero for 
each new case. 

= average total pressure over compressor face at 
time ”N". 

= mean value of FACT(N) for 50 previous values of 
N, when not zero. 

= standard deviation of FACT(N) for 50 previous 
values of N, when not zero. 



At the end of each 50 values of the index, a new mean and 
standard-deviation were computed and the test criterion was 
established for the next 50 values of the index. This simple 
test was used on all thirty-three runs made, and the results 
are summarized in the next section. Obviously, the test 
criterion could be lowered, with resulting greater occurrence 
of stall prediction, but this would also increase the number 
of false warnings. The value of test was initially set to 
predict stall for the index Case three without any false 

warnings. The test criterion could have been set to an opti- 
mum value for each particular index with possible better re- 
sults but this method was not used. The same test was applied 
to all cases. The statistical method has the real advantage 
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of not being subject to interpretation. All thirty-three 
computer plots are included in this thesis, since they cer- 
tainly are subject to personal interpretation as to the pre- 
diction of compressor stall. 

B. COMPARISON OF RESULTS 

There are two indicators of the results of the data re- 
duction. First are the computer plots, which provide an 
eyeball examination of results, A matrix summation of the 
computer plot results is provided by Table III, The second 
method was the use of statistical mean and standard-deviation 
discussed in the preceding section. A matrix summation of 
the statistical method is provided by Table IV. In both 
tables, where stall was predicted, the time to stall is given 
in milliseconds. 

There are several points of interest. As had been pre- 
viously expected, there is some correlation between the 
index and the Downs index, since both rely on circumferential 
distortion. Using the statistical method, all of the cases 
predicted by the index were also predicted by the Downs 

index. The Downs index predicted two additional cases uti- 
lizing this method. Using the computer plot interpretation, 
the method predicts six cases and the Downs index pre- 

dicts five, with four predictions common to both indices. 

If one were to use a combination of both vorticity based 
indices, the statistical method would have predicted eight 
out of the eleven cases presented. U.sing thi.s method there 
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TABLE III 

COMPUTER PLOT SUMT-URY 

(based on visual interpretation) 

RUN PRAT? & vraiTNSY SHOEMAKER DOVraS 

Kj,„ INDEX INDEX INDEX 


1 


No Prediction 


22.5 millisecs 


23.0 millisecs 


2 


3.0 millisecs 
3 same magnitude 


25.0 millisecs 


No Prediction 
one false peak 


3 


3.0 millisecs 


21 .0 millisecs 
4 false peaks 


6.0 millisecs 
2 false peaks 


4 


4.0 millisecs 
2 larger peaks 


No Prediction 


No Prediction 
one false peak 


5 


i-tU XxUUXUI/XUjJL 


3.0 rnillisccs 
2 smaller peaks 


No Prediction 


6 


No Prediction 


No Prediction 


No Prediction 


7 • 


No Prediction 


No Prediction 


No Prediction 


3 


4.5 millisecs 


No Prediction 


2.0 millisecs 


9 


4.0 millisecs 
one false peak 


l6.0 millisecs 
5 same magnitude 


2.0 millisecs 


10 


4.0 millisecs 


No Prediction 


2.0 millisecs 


1 1 


No Prediction 


6.0 millisecs 
2 larger peaks 


No Prediction 


totals 


6 Predictions 


6 Predictions 


5 Predictions 
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TABLE IV 

STATISTICAL Sin4>IARY 

RUN PRATT & VraiTNEY SHOEMAKER DOWNS 

K INDEX INDEX I^IDSX 


1 


No Prediction 


21.8 millisecs 
1 false warning 


22.28 millisecs 
2 false warnings 


2 


No Prediction 


23.9 millisecs 


No Prediction 
2 false warnings 


3 


0.9 millisecs 


20.2 millisecs 
1 false warning 


4,33 millisecs 
2 false warnings 


4 


No Prediction 


No Prediction 
1 false warning 


Mo Prediction 
1 false v^arning 


5 


No Prediction 


No Prediction 


No Prediction 


6 


No Prediction 


No Prediction 


No Prediction 


7 


No Prediction 


No Prediction 


19.33 millisecs 
1 false vrarning 


8 


2.88 millisecs 


No Prediction 
2 false warnings 


2.33 millisecs 


9 


0.4 millisecs 


35.3 millisecs 
2 false warnings 


0,0 millisecs 


10 


1 .62 millisecs 


No Prediction 
1 false warning 


2.33 millisecs 
1 false vrarning 


1 1 


No Prediction 


11.8 millisecs 


No Prediction 

_ _ i 


totals 


4 Predictions 


5 Predictions 


6 Predictions 
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would have been seventeen false warnings. The number of 
false warnings could be reduced by a reduction in the test 
criterion to a point where no loss in stall prediction would 
occur. 

There could be as much as seven milliseconds of time lag 
between actual stall time and predicted stall time. This 
presents a problem since the data could include the effects 
of the hammershock which occurs after stall, if the time lag 
is too great. It is certainly no problem to produce an index 
which will sense the hammershock and the question arises, did 
the index actually predict stall or were the effects of the 
hammershock being felt? By examining the numerical printout 
of average total pressure as a function of time and taking 
into account the possible seven-millisecond lag time, it was 
apparent that this could have been a problem in only four 
cases. These four cases should be held suspect and are cases 
3, 8, 9, and 10. It is noteworthy that the Shoemaker index 
predicted stall in two of these four cases sufficiently far 
in advance to preclude the problem. 
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VI. CONCLUSIONS 



An investigation into the problem of predicting compres- 
sor stall using instantaneous distortion parameters was made 
Some success was evidenced by both the statistical method 
and computer plot interpretation, but it was not possible to 
predict stall in every case using one or any combination of 
the indices explored. This seems to indicate that an axial 
compressor is subject to many instantaneous distortion pat- 
terns as well as to steady-state distortion, whereas one in- 
dex may only search for one category of distortion. If the 
cause of compressor stall was in every case instantaneous 
radial distortion over the outer portion of the compressor 
face, then in every case the Shoemaker index should predict 
stall. If the cause was in every case circumferential dis- 
trotion over some arbitrary sixty degree sector then the 
Downs index would be successful. If the distortion conforms 
to the ring pattern explored by the Pratt and Whitney index, 
then that method would succeed. In fact, it seems naive to 
believe that compressor stall could be predicted in every 
case using any one index. The indices explored in this the- 
sis and the indices used by industi’y today may certainly be 
improved upon, but it is doubtful that any one index will 
prove to be a panacea. 

In spite of the compound nature of compressor stall, in- 
stantaneous indices are useful. The Pratt and Whitney index 
used in this thesis was developed for the F-lllA with a 
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specific inlet-engine combination. It is in that context 
that distortion parameters would appear most useful - that 
is, to examine an inlet-engine system for compatibility. 

The performance of an inlet-engine system which too often 
experiences one type of distortion pattern may be improved 
upon; distortion indices would be very useful in isolating a 
predominant pattern. 

The use of vorticity as an index parameter has some real 
advantages. Both temperature and pressure gradients may be 
integrated into a single parameter. Vorticity may be visual- 
ized more simply and is a natural characteristic of the fluid. 
The Shoemaker index, when successful, harbingers stall suffi- 
ciently far in advance possibly to enable correction of the 
deficiency using engine controls. Vorticity appears to be a 
promising dimensionless parameter and deserves further inves- 
tigation . 
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APPENDIX A 



COMPUTER PLOTS 

The following thirty-three computer plots provide the 
instantaneous version of each of the three indices considered. 
In all cases, the horizontal scale is read in milliseconds. 

The vertical scale gives the index value and is a non-dimen- 
sional number. No consideration was given to index magnitude 
in predicting stall. For any given case, a strong peak just 
before stall in relation to that particular case was consid- 
ered to be sufficient for stall prediction. 
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Shoemaker Index 
case one 

index vs time in millisec 
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Dovais Index 
case one 

index vs time in millisec 



39 



I 



\ 



400 



FIGUHS A-3 



350 



300 



250 



200 



150 



ri: 



'I'M I' 



I ! 



' 



100 



iN'! 



' { 



IJil 



I* I iM!' 



1 II.; 



SlfiLL 



J: I. 

\ 



'lill 
I' M 



0.0 



50 



100 



150 



200 



250 



Pratt and Whitney Index 
case one 

index vs time in millisec 



40 



( \ 




Shoemaker Index 
case tv;o 

index vs time in millisec 
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Downs Index 
case two 

index vs tine in millisec 
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Pratt and Whitney Index 
case tvro 

index vs time in milliscc 
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Shoemaker Index 
case three 

index vs time in millisec 
' \ 
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Pratt and 'ATiitney Index 
case three 

index vs time in millisec 
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Shoemaker Index 
case four 

index vs time in millisec 
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Downs Index 
case four 

index vs tlrae in ailllsec 
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Pratt and V.Tiitney Index 
case four 

index vs time in inillisec 
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Shoemaker Index 
case five 

index vs time in raillisec 
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Downs Index 
case five 

index vs time In milliscc 



51 



' \ 




Pratt and vrhitney Index 
case five 

index vs time in raillisec 
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Shoemaker Index 
case six 

index vs time in millisec 
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Downs Index 
case six 

index vs time in millisec 
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Pratt and Whitney Index 
case six 



index vs time in raillisec 




Shoemaker Index 
case seven 

index vs time in raillisec 
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case seven 

. index vs time in millisec 
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Pratt and vrnltney Index 
case seven 

index vs time in millisec 
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Shoemaker Index 
case ci-;."ht 

index vs tine in millisec 
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Dovms Index 
case eight 

index vs time in inillisec 
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case eight 

index vs time in millisec 
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Shoemaker Index 
case nine 

index VO time in inillisec 
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Dovjns Index 
case nine 

index vs time in millisec 
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Pratt and >Thitney Index 
case nine 

index vs time in raillisec 
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Shoemaker Index 
case ten 

index vs time in inillisec 




Downs Index 
case ten 

index vs time In alllisec 
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Pratt and Whitney Index 
case ten 

index vs time in rnilliscc 
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case eleven 

index vs time in iiiillisec 
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Dovms Index 
case eleven 

index vs time in railllsec 
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Pratt and 'A'hitney Index 
case eleven 

index vs time in milliscc 
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SHUMAKER INDEX CASE 
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